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Abstract 

A new generation of neutrinoless double beta decay (Oz//3/3-decay) experiments with improved 
sensitivity is currently under design and construction. They will probe inverted hierarchy region 
of the neutrino mass pattern. There is also a revived interest to the resonant neutrinoless double- 
electron capture (Oi^ECEC), which has also a potential to probe lepton number conservation and 
to investigate the neutrino nature and mass scale. The primary concern are the nuclear matrix 
elements. Clearly, the accuracy of the determination of the effective Majorana neutrino mass from 
the measured O^/3/3-decay half-life is mainly determined by our knowledge of the nuclear matrix 
elements. We review recent progress achieved in the calculation of 0v{3f3 and Oi^ECEC nuclear ma- 
trix elements within the quasiparticle random phase approximation. A considered self-consistent 
approach allow to derive the pairing, residual interactions and the two-nucleon short-range correla- 
tions from the same modern realistic nucleon-nucleon potentials. The effect of nuclear deformation 
is taken into account. A possibility to evaluate Oi^/3/3-decay matrix elements phenomenologically is 
discussed. 

PACS numbers: 23.40.Hc, 21.60.Jz, 14.60.St, 12.60.Jv 

Keywords: 



1 



I. INTRODUCTION 



The physics community faces a challenging problem, finding whether neutrinos are indeed 
Majorana particles (i.e. identical to its own antiparticle) as many particle models suggest 
or Dirac particles (i.e. is different from its antiparticle). The best sensitivity on small 
Majorana neutrino masses can be reached in the investigation of neutrinoless double-beta 
decay (Oi//3/3-decay) [T|, , 

(A,Z) -> (A,Z + 2) + e- + e~ (1) 

and the resonant neutrinoless double-electron capture (Oz/ECEC) @, (sj, 

e ; + e; + (A,Z)^(A,Z-2y*. (2) 

A double asterisk in Eq. (J2J) means that, in general, the final atom (A, Z — 2) is excited with 
respect to both the electron shell, due to formation of two vacancies for the electrons, and 
the nucleus. Observing the Oz//3/3-decay and/or Oi/ECEC would tell us that the total lepton 
number is not a conserved quantity and that neutrinos are massive Majorana fermions. 

The search for the 0z//3/9-decay represents the new frontiers of neutrino physics, allowing 
in principle to fix the neutrino mass scale, the neutrino nature and possible CP violation 
effects. There are few tenths of nuclear systems Q, which offer an opportunity to study the 
0^/9/3-decay and the most favorable are those with a large Qpp-value. 

Neutrinoless double beta decay has not yet been found. The strongest limits on the half- 
life Ty 2 of the Oz//3/3-decay were set in Heidelberg-Moscow ( 76 Ge, 1.9 x 10 25 y) 0, NEM03 
( 100 Mo, 1.0 x 10 24 y) @, CUORICINO ( 130 Te, 3.0 x 10 24 y) @ and KamLAND-Zen ( 136 Xe, 
5.7 x 10 24 y) jsj experiments. However, there is an unconfirmed, but not refuted, claim 
of evidence for neutrinoless double decay in 76 Ge by some participants of the Heidelberg- 
Moscow collaboration]!] with half-life Ty 2 = 2.231 oJ 4 x 10 25 years. It is expected that the 



GERDA experiment [10J in a first phase will check this result relatively soon. 

The main aim of experiments on the search for 0i//9/9-decay is the measurement of the 
effective Majorana neutrino mass 

™W = H U lj m 3i ( 3 ) 

3 

where U e j is the element of Pontecorvo-Maki-Nakagawa-Sakata (PMNS) unitary mixing 
matrix and rrij is the mass of neutrino. For the most discussed case of mixing of three 
massive neutrinos (j=l,2,3) the PMNS matrix contains three charge parity (or CP) violating 
phases by assuming neutrinos to be Majorana particles. 

The effective Majorana neutrino mass can be calculated by using neutrino oscillation pa- 
rameters: an assumption about the mass of lightest neutrino, by chosing a type of spectrum 
(normal or inverted) and values of CP-violating phases. In future experiments, CUORE 
' 130 Te), EXO, KamLAND-Zen ( 136 Xe), MAJORANA ( 76 Ge), SuperNEMO ( 82 Se), SNO+ 
■'Nd), and others [H, H|, a sensitivity 



'150 ; 



\mpp\ ~ a few tens of meV (4) 

is planned to be reached. This is the region of the inverted hierarchy of neutrino masses. In 
the case of the normal mass hierarchy \mpp\ is too small, a few meV, to be probed in the 
O^/3/3-decay experiments of the next generation. 
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We note that it is reasonable to hope that the search for the Oz/ECEC of atoms, which are 
sufficiently long lived to conduct a practical experiment, may also establish the Majorana 
nature of neutrinos. This possibility is considered as alternative and complementary to 
searches for the Oz//3/3-decay. 

To interpret the data from the Oz//3/3-decay and the Oz/ECEC (neutrinoless double electron 
capture) accurately a better understanding of the nuclear structure effects important for the 
description of the nuclear matrix elements (NMEs) is needed. In this connection it is crucial 
to develop and advance theoretical methods capable to evaluate reliably NMEs, and to 
realistically assess their uncertainties. 



II. Oz//?/?-DECAY NMES: TWO-NUCLEON SHORT-RANGE CORRELATIONS 
AND UNCERTAINTIES 

The inverse value of the Oi//3/3-decay half-life for a given isotope (A, Z) can be written as 



1 rri/3/3 2 



1 l/2 



M'°T G°»{E ,Z). (5) 



Here, G 0u (E , Z) and M' 0l/ are, respectively, the known phase-space factor (E is the energy 
release) and the nuclear matrix element, which depends on the nuclear structure of the 
particular isotopes (A, Z), (A, Z + 1) and (A, Z + 2) under study. The phase space factor 
G 0u (E , Z) includes fourth power of unquenched axial- vector coupling constant g& and the 
inverse square of the nuclear radius R~ 2 , compensated by the factor R in M' 0l/ . The assumed 
value of the nuclear radius is R = tqA 1 ^ with r = 1.2 fm. 
The nuclear matrix element M' 0u is defined as 

M' 0u = M 0u . (6) 



Here, g e A is the quenched axial-vector coupling constant. This definition of M' 0u 11| allows 
to display the effects of uncertainties in g e /^ and to use the same phase factor G 0u (Eq, Z) 
when calculating the Oi//3/3-decay rate. 

The nuclear matrix element M 0u consists of the Fermi (F), Gamow- Teller (GT) and tensor 
(T) parts as 

M°» = -J^- + MX - 
(9A ff ) 2 



kl 



H F (r kl ) 

(7) 



(0, + | E^ + [-t4t^ + H GT{r k i)a kl - H T (r kl )S kl ]\0}) 



Here 



Ski = 3(<? fc ■ r kl )(ai ■ T k i) - o-fci, a k i = a k -ai. (8) 
The radial parts of the exchange potentials are 

H F,GT,T\f k i) = ~R / = dq. (9) 

it Jo q + E 
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where R is the nuclear radius and E is the average energy of the virtual intermediate 
states used in the closure approximation. The closure approximation is adopted in all the 
calculation of the NMEs relevant for Oz//3/3-decay with the exception of the QRPA. The 
functions h FGTT (q 2 ) are given by (l2| 



hGT{q ) - 3MO {g eff )2 4^ + 



P 

2 1 „4 



^ 3 q 2 + m 2 3 (g 2 + m 2 ) 2 / ' 

3 /1(g2) (V + "4) ~ (g 2 + m 2 ) 2 ) " 

For the vector normalized to unity and axial-vector form factors the usual dipole approx- 
imation is adopted: f v {q 2 ) = 1/(1 + q 2 /M 2 ) 2 , f A (q 2 ) = 1/(1 + q 2 /M\) 2 . My = 850 
MeV, and Ma = 1086 MeV. The difference in magnetic moments of proton and neutron 
is (/x p — fi n ) = 4.71, and qa = 1.254 is assumed. 

The above definition of the M® u includes contribution of the higher order terms of the nu- 
cleon current, and the Goldberger-Treiman PCAC relation, gp(q 2 ) = 2m p gA(q 2 )/(q 2 + m 2 ) 
was employed for the induced pseudoscalar term. 

The nuclear matrix elements M' 0u for the Ou/3/3 decay must be evaluated using tools of 
nuclear structure theory. Unfortunately, there are no observables that could be simply and 
directly linked to the magnitude of 0^/3/3 nuclear matrix elements and that could be used to 
determine them in an essentially model independent way. 

During many years two approaches were used: the Quasiparticle Random Phase Approx- 
imation (QRPA) dMi] and the Interacting Shell Model (ISM) [H] . There are substantial 



differences between both approaches. The QRPA treats a large single particle model space, 
but truncates heavily the included configurations. The ISM, by contrast, treats a small frac- 
tion of this model space, but allows the nucleons to correlate in many different ways. In the 
last few years several new approaches have been used for the calculation of the Oz//3/3-decay 
NMEs: the angular momentum Projected Hartree-Fock-Bogoliubov method (PHFB) [15], 
the Interacting Boson Model (IBM) |l6|, and the Energy Density Functional method (EDF) 
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The standard QRPA method consists of two steps. First, the mean field corresponding to 
the minimum of energy is determined and the like-particle pairing interaction is taken into 
account by employing the quasiparticle representation. In the second step the linearized 
equations of motion are solved in order to describe small amplitude vibrational-like modes 
around that minimum. In the renormalized version of QRPA (RQRPA) the violation of the 
Pauli exclusion principle is partially corrected. 

The drawback of QRPA is the fact that, unlike in BCS, the particle number is not 
conserved automatically, even on average. For realistic Hamiltonians the differences between 
averaged particle numbers on the RPA ground state and the exact particle numbers could be 
of the order of unity (an extra or missing neutron or proton). The selfconsitent renormalized 
QRPA method (SRQRPA) removes this drawback by treating the BCS and QRPA vacua 
simultaneously |l8|. For the neutron-proton systems the method was proposed and tested 
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TABLE I: Averaged Ou/3/3 nuclear matrix elements (M' 0u ) and their variance a (in parentheses) 
calculated within the QRPA and the RQRPA. Different types of two-nucleon short-range corre- 
lations (SRC) are considered: Milller-Spencer Jastrow SRCs (Jastrow) Fermi hypernetted 
chain SRCc (FHCh); unitary correlation operator method SRCs (UCOM) 121 ]; the coupled cluster 
method SRCs derived from the Argonne and CD-Bonn potentials 0- Threlsets of single particle 
level schemes are used, ranging in size from 9 to 23 orbits. The strength of the particle-particle 
interaction is adjusted so the experimental value of the 2^/3/3-decay nuclear matrix element is cor- 
rectly reproduced. Both free nucleon {g e / = 9A = 1.254) and quenched (ffl = 1 -0) values of 
axial- vector coupling constant are considered. We note that unlike in Refs. [111. E3. 25 ] tq = 1.2 
fm instead of ro = 1.1 fm is assumed. 



Nucleus 


eff 
9A 


meth. 






(M ) 






transition 








SRC 




CCM SRC 








Jastrow 


FHCh 


UCOM 


Argonne 


CD-Bonn 


lb Ge 


1.25 


QRPA 


4.92(0.19) 


5.15(0.17) 


5.98(0.27) 


6.34(0.29) 


6.89(0.35) 






RQRPA 


4.28(0.13) 


4.48(0.13) 


5.17(0.20) 


5.42(0.21) 


5.93(0.25) 




1.00 


QRPA 


4.18(0.15) 


4.36(0.15) 


4.97(0.23) 


5.20(0.22) 


5.63(0.27) 






RQRPA 


3.77(0.14) 


3.94(0.13) 


4.47(0.20) 


4.59(0.15) 


5.04(0.24) 


S2 Se 


1.25 


QRPA 


4.39(0.16) 


4.57(0.16) 


5.32(0.23) 


5.66(0.26) 


6.16(0.29) 






RQRPA 


3.81(0.14) 


3.97(0.14) 


4.59(0.17) 


4.84(0.21) 


5.30(0.22) 




1.00 


QRPA 


3.59(0.13) 


3.74(0.13) 


4.29(0.19) 


4.57(0.20) 


4.89(0.22) 






RQRPA 


3.17(0.10) 


3.32(0.10) 


3.79(0.13) 


4.00(0.15) 


4.29(0.16) 




1.25 


QRPA 


1.22(0.03) 


1.23(0.04) 


1.77(0.02) 


2.07(0.10) 


2.28(0.03) 






RQRPA 


1.31(0.15) 


1.33(0.15) 


1.77(0.02) 


2.01(0.17) 


2.19(0.22) 




1.00 


QRPA 


1.32(0.08) 


1.34(0.07) 


1.73(0.10) 


1.90(0.12) 


2.11(0.12) 






RQRPA 


1.22(0.12) 


1.24(0.12) 


1.57(0.14) 


1.69(0.13) 


1.88(0.16) 


inn i r 


1.25 


QRPA 


3.64(0.21) 


3.73(0.21) 


4.71(0.28) 


5.18(0.36) 


5.73(0.34) 






RQRPA 


3.03(0.21) 


3.12(0.21) 


3.88(0.26) 


4.20(0.34) 


4.67(0.31) 




1.00 


QRPA 


2.96(0.15) 


3.02(0.15) 


3.74(0.21) 


4.03(0.27) 


4.44(0.24) 






RQRPA 


2.55(0.13) 


2.63(0.13) 


3.20(0.17) 


3.43(0.25) 


3.75(0.21) 




1 9^ 






Q 1 1 (C\ 01 "| 


q j Ala 1 9"i 


3.86(0.29) 


4.35(0.16) 






RQRPA 


2.64(0.17) 


2.75(0.19) 


3.21(0.22) 


3.34(0.24) 


3.72(0.26) 




1.00 


QRPA 


2.38(0.17) 


2.47(0.17) 


2.88(0.17) 


2.99(0.23) 


3.31(0.21) 






RQRPA 


2.14(0.14) 


2.21(0.14) 


2.55(0.17) 


2.69(0.19) 


2.92(0.21) 


128 Te 


1.25 


QRPA 


3.97(0.14) 


4.15(0.15) 


5.04(0.15) 


5.38(0.17) 


5.99(0.17) 






RQRPA 


3.52(0.13) 


3.68(0.14) 


4.45(0.15) 


4.71(0.17) 


5.26(0.16) 




1.00 


QRPA 


3.11(0.09) 


3.23(0.10) 


3.88(0.11) 


4.11(0.13) 


4.54(0.13) 






RQRPA 


2.77(0.09) 


2.88(0.09) 


3.44(0.10) 


3.62(0.12) 


4.00(0.12) 


130 Te 


1.25 


QRPA 


3.56(0.13) 


3.72(0.14) 


4.53(0.12) 


4.77(0.15) 


5.37(0.13) 






RQRPA 


3.22(0.13) 


3.36(0.15) 


4.07(0.13) 


4.27(0.15) 


4.80(0.14) 




1.00 


QRPA 


2.55(0.08) 


2.93(0.08) 


3.52(0.07) 


3.69(0.11) 


4.11(0.08) 






RQRPA 


2.15(0.14) 


2.66(0.09) 


3.17(0.08) 


3.29(0.11) 


3.69(0.09) 


136 Xe 


1.25 


QRPA 


2.16(0.13) 


2.25(0.12) 


2.73(0.13) 


2.88(0.14) 


3.23(0.14) 






RQRPA 


2.02(0.12) 


2.11(0.14) 


2.54(0.15) 


2.68(0.16) 


3.00(0.17) 




1.00 


QRPA 


1.70(0.09) 


1.77(0.09) 


2.12(0.11) 


2.21(0.10) 


2.47(0.09) 






RQRPA 


1.59(0.09) 


1.66(0.10) 


1.97(0.11) 


2.06(0.11) 


2.30(0.12) 
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on the exactly solvable simplified models in Ref . 19[ . It is a generalization of the procedure 
proposed earlier in (20| . 

In the QRPAthe phonon operators are defined as 

Q^J,M — ^pn[-^^m)J^\pn)J,M ~ ^{pn) jA(pn) J,m] , (H) 

where X^^j and Y^ pn )j are the usual variational amplitudes, and ^4( pn ) JM is the angular 
momentum coupled two-quasiparticle creation operator, p, n signify the quantum numbers 
of the proton, respectively neutron, orbits. The X and Y amplitudes, as well as the corre- 
sponding energy eigenvalues Uk are determined by solving the QRPA eigenvalue equations 
for each J n 

(-\- B a) (?)-(?)■ (12) 

The matrices A and B above are determined by the Hamiltonian rewritten in terms of the 
quasiparticle operators: 

Ktfn> = (O|(aH) (JM)t ^(4'40 (JM) |O) (13) 
B J m , p , n , = (0|^(atat) (J - M) (-l) M («>l) (JM) |0) 

Here, \0) is the BCS vacuum state. 

In the RQRPA and SRQRPA instead of \0) the correlated QRPA ground state |0q RPA ) 
is considered. Then instead of the standard X and Y everywhere and also in the QRPA 
equations of motion the renormalized amplitudes are used: 

'Y m — T) 1 / 2 Y m V" 1 — T> 1 / 2 V m ('\A\ 

yv {pn,J 7 ') ~ ^pn yv (pn,J^)i 1 {pn,.J 7T ) ~ ^pn 1 (pn,J 7r )5 \ L ^J 

where renormalization factors T> pn are given by 

T^pn = (^QRPA\[^-(pn)J,M^ ^(pn) J,m] I OqHPa) = 1 ~ £p ~ £n (15) 

= 1 - ^^^V pn \^ k {2J + l)\Y J p t\ 2 ) 

- ^T-T^p'Vp'n (Ej,*(2J + i)\Y J / n \ 2 ) ■ 
Here, £ n ( p ) is the expectation value of the number of quasiparticles in the orbit n(p), 

(®QRPa\ a n(p) a n(p) \0qRPa) 

Zn(p) = , = • (16) 

\j l 3n{p) + 1 

a ta( P ),mi a jn( P ),m are the creation and annihilitation operators for the quasiparticle with quan- 
tum numbers n(p),m. The renormalizattion coefficients V pn and the quasiparticle occupa- 
tion numbers £j can be obtained iteratively using the equations of motion of the (S)RQRPA. 

In the correlated QRPA ground state the occupation numbers are no longer the pure 
BCS quantities. Instead, they depend, in addition, on the solutions of the QRPA equations 
of motion for all multipoles J, and can be evaluated using 

n n(p) = (QQRPA\^mCn(p),m C n(p),m\®QRPA) (17) 
~ {2j n(p) + 1) [4) + - ujL)) Cn(p) 



6 



Here, c- m is the creation operator for a proton in the orbit j p or a neutron in the orbit j n and 
Cj t m is the corresponding annihilation operator. The amplitudes Vj p and Vj n are obtained by 
solving the gap equations. 

In SRQRPA the BCS equations are reformulated. This is achieved by recalculating 
the u and v amplitudes from the minimum condition of the RQRPA ground-state energy. 
In SRQRPA thus the state around which the vibrational modes occur is no longer the 
quasiparticle vacuum, but instead the Bogoliubov transformation is chosen is such a way 
that provides the optimal and consistent basis while preserving the form of the phonon 
operator, Eq. ffTTT) . 

In practice, the SRQRPA equations are solved double iteratively. One begins with the 
standard BCS u, v amplitudes, solves the RQRPA equations of motion and calculates the 
factors D pn . The u, v amplitudes are recalculated and the procedure is repeated until the 
selfconsistency is achieved. Numerically, the double iteration procedure represents a chal- 
lenging problem. It was resolved in [2lj where instead of the G-matrix based interaction the 
pairing part (and only that part) of the problem was replaced by a pairing interaction that 
uses a constant matrix element whose value was adjusted to reproduce the experimental 
odd-even mass differences. 

In the QRPA, RQRPA, and SRQRPA the M 0u is written as the sum over the virtual 
intermediate states, labeled by their angular momentum and parity J n and indices ki and 
k f : 

m k = E (-iy n+j p' +J+J V2j+r{ h 3n I}* 



{p(l),p'(^J || f(r 12 )0 K f(r 12 ) || n(l),n'(2);J) x 
(0;il[45„']j||J 7r A; / )(J-A; / |J^)(J 7r fc i ||[ C +g n ] J ||0+) . 



(18) 



The operators Ok,K = Fermi (F), Gamow- Teller (GT), and Tensor (T), contain neutrino 

potentials and spin and isospin operators, and RPA energies Ej" 1 . Two separate multipole 
decompositions are built into Eq. (|T8|) . One is in terms the J n of the virtual states in the 
intermediate nucleus, the good quantum numbers of the QRPA and RQRPA. The other 
decomposition is based on the angular momenta and parities J 1 * of the pairs of neutrons 
that are transformed into protons with the same J 71 ' . The nucleon orbits are labeled in 
Eq. (TTHT) by p, p', n, n'. 

The QRPA-like approaches do not allow to introduce short-range correlations (SRCs) 
into the two-nucleon relative wave function. The traditional way is to introduce an ex- 
plicit Jastrow-type correlation function f(r 12) into the involved two-body transition matrix 
elements (see Eq. [TBI . In the parametrization of Miller and Spencer [22] we have 

/(r 12 ) = 1 - ce~ ar \l - br 2 ), a = 1.1 fni' 2 , b = 0.68 frnT 2 . (19) 

These two parameters (a and b) are correlated and chosen in the way that the norm of the 
relative wave function is conserved. 

Recently, it was proposed [l3j to adopt instead of the Jastrow method the UCOM ap- 
proach for description of the two-body correlated wave function j23[ . The UCOM method 
produces good results for the binding energies of nuclei already at the Hartree-Fock level 
2!. 
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TABLE II: The calculated ranges of the nuclear matrix element M 0v evaluated within the QRPA 
(column 2), RQRPA (column 4) and SRQRPA (column 6), with standard (g e J f = g A = 1.254) 
and quenched (g A = 1.0) axial- vector couplings and with the coupled cluster method (CCM) 
CD-Bonn and Argonne short-range correlation (SRC) functions. Columns 3, 5 and 7 give the 
Oi^/3/3-decay half-life ranges corresponding to values of the matrix-elements in columns 2, 4 and 6 
for = 50 meV. T®J 2 ~ exp is the experimental lower bound on the Oz^/3/3-decay half-life for a 

given isotope. 



Nucl. QRPA 


RQRPA 


SRQRPA 




M°» Ty 2 [y] 


M'»» T»J 2 [y] 


M'"» Ty 2 [y] 





Ge (5.0,7.2) (3.0,6.3) x 10 26 (4.5,6.2) (4.1,7.9) x 10 26 (4.3,6.2) (4.0,8.6) x 10 26 1.9 x 10 25 [5] 

Se (4.4,6.4) (8.5,18.) x 10 25 (3.8,5.6) (1.2,2.4) x 10 26 (3.9,6.1) (9.5,22.) x 10 25 3.2 x 10 23 [6] 

°Mo (3.7,6.1) (5.9,15.) x 10 25 (3.2,5.0) (8.8,21.) x 10 25 (4.0,5.5) (7.3,13.) x 10 25 1.0 x 10 24 [6] 

°Te (3.6,5.5) (7.4,18.) x 10 25 (3.2,4.7) (1.0,2.2) x 10 26 (3.6,5.1) (8.5,17.) x 10 25 3.0 x 10 24 [7] 

6 Xe (2.1,3.4) (1.9,4.8) x 10 26 (2.0,3.2) (2.1,5.5) x 10 26 (2.4,3.6) (1.6,3.7) x 10 26 5.7 x 10 24 [8] 



A self-consistent calculation of the 0i//9/9-decay NMEs in the QRPA-like approaches was 



performed in [25]. The pairing and residual interactions as well as the two-nucleon short- 
range correlations were for the first time derived from the same modern realistic nucleon- 
nucleon potentials, namely from the charge-dependent Bonn potential (CD-Bonn) and the 
Argonne V18 potential. A method of choice was the coupled cluster method (CCM) [2fj|. For 
purpose of numerical calculation of the 0z//3/3-decay NMEs the CCM short-range correlation 
functions were presented in an analytic form of Jastrow-like function as [25] 

/A,B(r 12 ) = l - ce- ar \l-br 2 ). (20) 

The set of parameters for Argonne and CD-Bonn NN interactions is given by 

f A (r 12 ) : a = 1.59 fm~ 2 , b = 1.45 fm~ 2 , c = 0.92, 
f B (r 12 ) : a = 1.52 fm~ 2 , b = 1.88 fm~ 2 , c = 0.46. 

(21) 

The calculated NMEs with these short-range correlation functions agree within a few per- 
centages with those obtained without this approximation. We note that the dependence of 
the SRC on the value of oscillator length b is rather weak. 

In Table fl] the QRPA and RQRPA results are presented separately for different types of 
two-nucleon short-range correlations (SRC) are considered: Milller- Spencer Jastrow SRCs 
(Jastrow) [lH; Fermi hypernetted chain SRCc (FHCh); unitary correlation operator method 
SRCs (UCOM) Il2j; the coupled cluster method SRCs derived from the Argonne and CD- 
Bonn potentials [25J based on an extension of the Brueckner theory (Coupled Cluste Method 
= CCM). Two different values of the axial coupling constant, free nucleon gf f = g A = 1.254 
and quenched g A = 1.0, are taken into account. The strength of the particle-particle 
interaction is adjusted so the experimental value of the 2i//9/3-decay nuclear matrix element 
is correctly reproduced [ll|]. The NME calculated within this procedure, which includes 
three different model spaces, is denoted as the averag ed 0^,5-decay NME (M ,0u ). We note 
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that the values of NMEs become essentially independent on the size of the single-particle 
basis and rather stable with respect to the possible quenching of the qa- 

From Table [T]it follows that the QRPA values are about 10-15% larger in comparison with 
the RQRPA values. The largest NMEs are those calculated with the CCM CD-Bonn corre- 
lation function. In comparison with them the NMEs obtained with the CCM CD-Argonne 
correlation function and the UCOM SRCs are about 10% smaller. This is explained by the 
fact that the CCM Argonne correlation function cuts out more the small r 12 part from the 
relative wave function of the two-nucleons as the CCM CD-Bonn correlation function. The 
smallest in magnitude are matrix elements for the Ovflfl decay obtained with the traditional 
approach of using the Miller-Spencer Jastrow SRC and the Fermi hypernetted chain SRCc. 

In Table [Til we show the calculated ranges of the nuclear matrix element M 0u evaluated 
within the QRPA, RQRPA Q and SRQRPA Q in a self-consistent way with the CCM 
CD-Bonn and Argonne SRC functions by assuming both the standard (g^ = 1.254) and 
quenched (gA = 1-0) axial- vector couplings. These ranges quantify the uncertainty in the 
calculated 0i//9/3-decay NMEs of a given QRPA-like approach. By comparing the SRQRPA 
with the RQRPA results we conclude that the requirement of conserving the particle number 
have not caused substantial changes in the value of the 0vf3f3 matrix elements in that case. 

Given the interest in the subject, in Table HT1 we show also the range of predicted 0^/3(3- 
decay half-lives of 76 Ge, 82 Se, 100 Mo, 130 Te and 136 Xe corresponding to full range of M' 0u for 
\ m pp\ = 50 meV. This is a rather conservative range within the considered QRPA framework. 
It represents roughly a required sensitivity of the O^/3/3-decay experiment in the case of 
inverted hierarchy of neutrino masses, which can be compared with the current bound on 
the Oi//3/3-decay half-life Ty^~ exp . 



III. /3/3 DECAY OF DEFORMED NUCLEI WITHIN QRPA 



One of the best candidates for searching 0^/3/3 decay is 150 Nd since it has the second 
highest endpoint, Qpp =3.37 MeV, and the largest phase space factor for the decay (about 
33 times larger than that for 76 Ge, see e.g. |28|). The SNO+ experiment at the Sudbury 
Neutrino Observatory will use a Nd-loaded scintillator to search for neutrinoless double beta 
decay by looking for a distortion in the energy spectrum of decays at the endpoint [29]. 

However, 150 Nd is well known to be a rather strongly deformed nucleus. This strongly 
hinders a reliable theoretical evaluation of the corresponding Oi/S/3-decay NMEs (for in- 
stance, it does not seem feasible in the near future to reliably treat this nucleus within the 
large-scale nuclear shell model (ISM), see, e. g., R ef. (13] )• Recently, more phenomenolog- 
ical approaches like the pseudo-SU(3) model [30], the PHFB approach [151 ]. the IBM [16] , 
and the EDF [l7[ have been employed to calculate M 0u for strongly deformed heavy nuclei 
(a comparative analysis of different approximations involved in some of the models can be 
found in Ref. (HI). The results of these models generally reveal a substantial suppression of 
M 0u for i50 Nd ag comparec i wrt h the QRPA result of Ref. [Tl| where 150 Nd and 150 S m were 
treated as spherical nuclei. However, the calculated NMEs M 0u for 150 Nd show a rather 
significant spread. 

One of the most up-to-date microscopic ways to describe the effect of nuclear deforma- 
tion on /3/3-decay NMEs M 2u and M 0u is provided by the QRPA. A QRPA approach for 
calculating /3/3-decay amplitudes in deformed nuclei has been developed in a series of pa- 



pers 



32 



in Ref. 



3Q. M 2u were calculated in Refs. j32|, with schematic separable forces, and 
with realistic residual interaction. It was demonstrated in Refs. [32-34 that 



34 
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deformation introduces a mechanism of suppression of the M 2u matrix element which gets 
stronger when deformations of the initial and final nuclei differ from each other. A similar 
dependence of the suppression of both M 2v and M 0u matrix elements on the difference in 
deformations has been found in the PHFB 
In Refs. 
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35. 36 



the first QRPA calculations of M 0l/ with an account for nuclear defor- 
mation were done. The calculations showed a suppression of M 0u for 150 Nd by about 40% 
as compared with our pre vious spherical QRPA result I n the next section we review 
the results of Refs. 32^36 . 



A. Formalism 

The NMEs M 2v and M 0u , as the scalar measures of the decay rates, can be calculated in 
any coordinate system. For strongly deformed, axially symmetric, nuclei the most convenient 
choice is the intrinsic coordinate system associated with the rotating nucleus. This employs 
the adiabatic Bohr-Mottelson approximation that is well justified for 150 Nd, 160 Gd and 160 Dy, 
which indeed reveal strong deformations. As for 150 Sm, the enhanced quadrupole moment 
of this nucleus is an indication for its static deformation. 

Though it is difficult to evaluate the effects beyond the adiabatic approximation employed 
here, one might anticipate already without calculations that the smaller the deformation is, 
the smaller should be the deviation of the calculated observables from the ones obtained in 
the spherical limit. In this connection it is worth noting that spherical QRPA results can 
exactly be reproduced in the present QRPA calculation by letting the deformation vanish, 
in spite of the formal inapplicability of the adiabatic ansatz for the wave function in this 
limit. 

We give here for completeness the formalism of the QRPA calculations of NMEs M 2v 
and M 0u in deformed nuclei as developed in Refs. 32 36|. 



Nuclear excitations in the intrinsic system \K n ) are characterized by the projection of 
the total angular momentum onto the nuclear symmetry axis K (the only projection which 
is conserved in strongly deformed nuclei) and the parity ir. 

The intrinsic states \K*,m) are generated within the QRPA by a phonon creation oper- 
ator acting on the ground-state wave function: 

iK^jm) = Q mt K\Qg.s.)'i Qm,K = ^jm. K^pn. K ~ Y™ K A pn ^ K . (22) 

pn 

Here, A pnK = a p a n t K = apan are the two-quasiparticle creation and annihilation 

operators, respectively, with the bar denoting the time-reversal operation. The quasiparticle 
pairs pn are defined by the selection rules Q p — Q n = K and ir p 7T n = n, where tt t is the single- 
particle (s.p.) parity and Q T is the projection of the total s.p. angular momentum on the 
nuclear symmetry axis (r = p,n). The s.p. states \p) and |n) of protons and neutrons are 
calculated by sol ving the Schrodinger equation with the deformed axially symmetric Woods- 
Saxon potential [34J. In the cylindrical coordinates the deformed Woods-Saxon s.p. wave 
functions |rO T ) with Q T > are decomposed over the deformed harmonic oscillator s.p. 
wave functions (with the principal quantum numbers (Nn z A)) and the spin wave functions 
|S = ±i): 

|<rfi T ) = £ b Nnz x\Nn z A r = n T -Y;}\Y;}, (23) 

JVn.S 
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where N = n± + n z (n± = 2n p + |A|), n z and n p are the number of nodes of the basis 
functions in the z- and p-directions, respectively; A = Q — X and £ are the projections of 
the orbital and spin angular momentum onto the symmetry axis z. For the s.p. states with 
the negative projection Q T = —\Q T \, which are degenerate in energy with Q T = |0 T |, the 
time- reversed version of Eq. (I2"3"j) is used as a definition (see also Ref. [34]]). The states (r, f ) 
comprise the whole single-particle model space. 

The deformed harmonic oscillator wave functions \Nn z A) can be further decomposed 
over the spherical harmonic oscillator ones \n r lA) by calculating the corresponding spatial 
overlap integrals A = (n r lA\Nn z A) (n r is the radial quantum number, I and A are the 
orbital angular momentum and its projection onto z-a.xes, respectively), see Appendix of 
Ref. [34j for more details. Thereby, the wave function fl23|) can be reexpressed as 

\rQ T ) = £B;|?^ T >, (24) 

V 

where \rjfl T ) = Z)C/q T _ e 1 ^\n r lA = Q T — E)|E) is the spherical harmonic oscillator wave 

S T 2 

function in the j-coupled scheme (77 = (n r lj)), and = I]C / - ? ^ r _ s 1 „ A^ x q t _ s bjsr nz x, 

S T 2 

with C° { n T _ E 1 E being the Clebsch-Gordan coefficient. 
The QRPA equations: 



A{K) B(K) \ (X$\ (X% } 
-B(K) -A(K) J \Yg) ~ UK ' m { J ' 

with realistic residual interaction are solved to get the forward X^, backward amplitudes 
and the excitation energies uf^ and uf^ of the m-th K w state in the intermediate nucleus. 
The matrix A and B are defined by 

Apn )P ' n '(K) = 6 pn ^i n i(E p + E n ) + g P p(u p u n Up'U n > + v p v n Vpiv n i)Vp^w -1 



- 9 P h(u p v n Up'V n ' + v p u n Vp>u n >)Vp„>..j 



pnp'n' 
pn'p'n 



B P n t p' n '(K) — —gpp{upU n Vpiv n i + VpV n Up'U n i)V p ^ r ,' 



pnp'n' 

-gph{u P v n Vp>v n , + VpU n u pl v nl )Vp nl pi n (26) 

where E p + E n are the two-quasiparticle excitation energies, V pntP > n ' and V pnp ,-, are the p — h 
and p — p matrix elements of the residual nucleon-nucleon interaction V, respectively, u T 
and v T are the coefficients of the Bogoliubov transformation. 

As a residual two-body interaction we use the nuclear Brueckner G matrix, which is a 
solution of the Bethe-Goldstone equation, derived from the charge-depending Bonn (Bonn- 



CD) one boson exchange potential, as used also in the spherical calculations of Ref. [TJj . The 
G matrix elements are originally calculated with respect to a spherical harmonic oscillator 
s.p. basis. By using the decomposition of the deformed s.p. wave function in Eq. ([Ml), the 
two-body deformed wave function can be represented as: 

\pn) = E F p f pn J V p Vn ,JK), (27) 

where \r} p r} n ,JK) = Y. mvm Si% v u mn Vlv m v)\ r lnm n ) , and E^ nr]n 
BP p B™ n (-l) jn ~ nn Cj p Q pjn ^ nn is defined for the sake of simplicity ((-l)J"~ n ™ i s the phase 
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arising from the time-reversed states \n)). The particle-particle V pfl ^ p ,~i and particle-hole 
V pn r ) p'n interaction matrix elements in the representation (|26|) for the QRPA matrices A, B 
( 125|) in the deformed Woods-Saxon single-particle basis can then be given in terms of the 
spherical G matrix elements as follows: 

V pn, p'n' =~ 2 EE E F p^nr, n F p>n p m'n n , G ( r lp r lnrip'riri,J), (28) 

J VpVnVp'Vn' 

JK' JK' 

Vpn', p'n = 2^^^^^^"^'^ ( 29 ) 

J IpTlnVp'Vn' 

where K' pn , = Q p + Q n > = Q p > + Q„. 

The structure of the intermediate |0 + ) and |1 + ) states is only needed within the QRPA 
to calculate 2i//3/3-decay NMEs M 2u [34], whereas all possible \K*) states are needed to 
construct the NMEs M 0u . 

The matrix element M^ T is given within the QRPA in the intrinsic system by the following 
expression: 

M 2 J T E E < /l^l K+ ' m ^^ + ' m /l ir+ ' m ^^ + ' m -l^' ^ (go) 

Instead of the usual approximation of the energy denominator in Eq. ( 1301) as CjK,m imf = 
{^K,m s + ^,mJ/2 (see, e.g., Refs. |32|, |33|), here another prescription is used in which 
the whole calculated QRPA energy spectrum is shifted in such a way as to have the first 
calculated 1 + state exactly at the corresponding experimental energy. In this case the energy 
denominator in Eq. (J30]) acquires the form ujK, mi m f = (uK,m f —UK,i J +UK,m i — UK,u)l'2>+u 1 +, 
with Q x + being the experimental excitation energy of the first 1 + state measured from the 
mean g.s. energy (E 0i + E 0f )/2. 

The two sets of intermediate nuclear states generated from the initial and final g.s. do 
not come out identical within the QRPA. Therefore, the overlap factor of these states is 
introduced in Eq. (!30|) [32] as follows: 

(K\ m f \K\ mi ) = E [K/kKk - Y$Y&) {BCS f \BCSi). (31) 

ki f 

The factor IZi^, which includes the overlaps of single particle wave functions of the initial 
and final nuclei is given by: 

n w = (pp p \p' Ppl ){ufuf + v^PXnpnln'p^Ku^ + v$>v%>), (32) 

and the last term (BCSf\BCSi) in Eq. fl3Tl) corresponds to the overlap factor of the initial 
and final BCS vacua in the form given in Ref. [3~2| . 

The matrix element M 0u is given within the QRPA in the intrinsic system by a sum of 
the partial amplitudes of transitions via all the intermediate states K 71 : 

M 0u = E M 0u {K n ) ; M 0v (K*) = E s^ ef) O a (K*). (33) 

Here, we use the notation of Appendix B in Ref. 12j], a stands for the set of four single- 
particle indices {p,p',n,n'}, and O a (K n ) is a two-nucleon transition amplitude via the K n 
states in the intrinsic frame: 

O a (K 7r ) = E (0;|cJc n |^m / )(^m / |^m J )(^m i |4,c n ,|0+). (34) 

mi,m f 
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The two sets of intermediate nuclear states generated from the initial and final g.s. (la- 
beled by rrii and mf, respectively) do not come out identical within the QRPA. A standard 
way to tackle this problem is to introduce in Eq. (|34|) the overlap factor of these states 
(K 71 'm f\K^ 'mi) , whose representation is given below, Eq. ( 13 7p . Two-body matrix elements 
s (def) Q f neu t r i no potential in Eq. (I3"3"|) in a deformed Woods-Saxon single-particle basis 
are decomposed over the the spherical harmonic oscillator ones according to Eqs. (j27|29p : 

(def) \ - \p F JK F JK Jsph) I j\ 

°pp'nn' Z_^i PVpnrfn p t ri p in l 'q n i J ri v 'q p iri n ri n i\' J )■> \ UKJ J 

3 VpVp> 

VnVn' 



( s P h ) 
^pp 



tAJ) = E(-l) J " +V+J+ ^{^, f p , ^}(p(l),p , (2); l 7||a(l,2)||n(l),n'(2); l 7)(36) 



where J = \/2J + 1, and Ot(l, 2) is the neutrino potential as a function of coordinates of 
two particles, with £ labeling its Fermi (F), Gamow- Teller (GT) and Tensor (T) parts. 

The particle-hole transition amplitudes in Eq. ( )3^|) can be represented in terms of the 
QRPA forward and backward amplitudes along with the coefficients of the Bogoli- 
ubov transformation u T and v T [34| : 

(Oj |CpC n |-K" mf) = V p U n X p ^ K TT + UpDnYpn^-n , 
(K mi |CpC n |0^) = UpV n X pl l K TT + VpUnYpn'xn. 

The overlap factor in Eq. can be written as: 

(K*m f \K*m t ) = Y\ X ? s K«Kk* ~ Y^Y^}TZ lfh {BCS f \BCS,) (37) 
klf 

Representations for TZy^ and the overlap factor (BCSf\BCSi) between the initial and final 
BCS vacua are given in Ref. 32| . 



B. Calculation results 

The NMEs M 2v and M 0l/ were calculated according to the above formalism in Refs. (3 



361 ] . These articles contain detailed description of the choice of the model parameters and 
comparison between different approximations. Here we only briefly repeat the key points of 
the calculations. 

Only quadrupole deformation is taken into account in the calculations [32-361. The single- 



particle Schrodinger equation with the Hamiltonian of a deformed Woods-Saxon mean field 
is solved on the basis of an axially-deformed harmonic oscillator. Decomposition of the 
obtained deformed single-particle wave functions is performed over the spherical harmonic 
oscillator states within the seven lowest major shells. The geometrical quadrupole defor- 
mation parameter 02 of the deformed Woods-Saxon mean field is obtained by fitting the 
experimental deformation parameter (3 = \J\%p-> where r c is the charge rms radius and Q p 
is the empirical intrinsic quadrupole moment. The experimental values of j3 can be derived 
from the laboratory quadrupole moments measured by the Coulomb excitation reorienta- 
tion technique, or from the corresponding B(E2) values [39[. Experimental values extracted 
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TABLE III: Calculated NMEs M' 0l/ and corresponding decay half-lives (assuming m«g=50 meV) 

721 ] within the 



for Oi//3/3-decays of 76 Ge, ibU Nd, iBU Gd 



150t 



160/ 



36j, and for Oi/ECEC of 152 Gd, 164 Gd, 180 W 
deformed QRPA. The results obtained for realistic deformations of the nuclei are labeled as "def , 
whereas those obtained in the spherical limit, i.e. fa = 0, are labeled as "sph". 



Nuclear transition 


eff 
9a 




sph (p 2 =0J 




ciei 










r i/2 Ly r J 


M 


r i/ 2 iy r \ 










\J is LJ LJ llvv. U y 








76 Ge^ 76 Se 


0.94 


4.10 


9.4 x 10 26 


4.00 


9.8 x 10 26 






1.25 


5.30 


5.6 x 10 26 


4.69 


7.2 x 10 26 




150 Nd^ 150 Sm 


0.94 


4.52 


2.3 x 10 25 


2.55 


7.1 x 10 25 






1.25 


6.12 


1.2 x 10 25 


3.34 


4.1 x 10 25 




160 Gd ^160 Dy 


0.94 




Oi/ECEC 


3.76 


2.3 x 10 26 




152 Gd^ 152 Se (KLi) 


1.269 


7.50 


(8.7 x 10 27 , 8.9 x 10 29 ) 


3.23 


(4.7 x 10 28 , 4.8 x 


10 29 ) 


i64 Er ^i64 Dy ( LiLi ) 


1.269 


7.20 


(1.0 x 10 32 , 1.1 x 10 32 ) 


2.64 


(7.5 x 10 32 , 8.4 x 


10 32 ) 


180 W ^180 Hf (KK) 


1.269 


6.22 


(1.4 x 10 30 , 2.0 x 10 30 ) 


2.05 


(1.3 x 10 31 , 1.8 x 


10 31 ) 



from the B(E2) have smaller experimental eroors. But deformations extracted from the re- 
orientation effect are in principle the better values, but have large errors. The fitted values 
of the parameter /3 2 of the deformed Woods-Saxon mean field, which allow us to reproduce 
the experimental /3, are listed in Table 1 of Ref. j36[. The spherical limit, i.e. 02 = 0, is 
considered as well, to compare with the earlier results of Ref. [Hj]. 

The nuclear Brueckner G matrix, obtained by a solution of the Bethe-Goldstone equa- 
tion with the Bonn-CD one boson exchange nucleon-nucleon potential, is used as a residual 



two-body interaction in Refs. [34H36| . Then the BCS equations are solved to obtain the 
Bogoliubov coefficients, gap parameter and chemical potentials. To solve the QRPA equa- 
tions, one has to fix the particle-hole g p h and particle-particle g pp renormalization factors 
of the residual interaction, Eqs. (I26j) . A value of g p h = 0.90 was determined by fitting the 
experimental position of the Gamow- Teller giant resonance (GTR) for 76 Ge. The same value 
of g p h wa s then used for all nuclei in question, and led to a good fit to the experimental 
GTR energy for 150 Nd, measured very recently (37| . 

The parameter g pp can be determined by fitting the experimental value of the 2vf3f3- 
decay NMEs M"q T MeV -1 [38| for each nucleus in question. To account for the quenching 
of the axial- vector coupling constant gA, the quenched value g e j[ = 0.75 -g& was used in the 
calculation along with the bare value g& = 1.25. The quenching factor of 0.75 comes from a 
recent experimental measurement of GT strength distribution in 150 Nd [37| . The two sets of 
the fitted values of q p „ corresponding to the cases without or with quenching of g\ are listed 
in Table 1 of Ref. [3Gf] . Note, that this fitting procedure leads to realistic values g pp ~ I. 

Having solved the QRPA equations, the two-nucleon transition amplitudes (1541) are calcu- 
lated and, by combining them with the two-body matrix elements of the neutrino potential, 
the total Oz//3/3 NMEs M Qv (1331) is formed. Such a computation is rather time consuming 
since numerous programming loops are needed to calculate the decompositions of the two- 
body matrix elements in the deformed basis over the spherical ones. Therefore, to speed up 
the calculations the mean energy of 7 MeV of the intermediate nuclear excitation energies is 
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TABLE IV: The matrix elements M 0u for the 0v(3(3 decay 150 Nd-» 150 Sm calculated in different 

1/2 



models. The corresponding half-lives Ty 2 ( m years) for an assumed effective Majorana neutrino 



mass (mpp) = 50 meV are also shown. 



Method def. QRPA [36] pseudo-SU(3) [30] PHFB [15] IBM [16] EDF [17] 



M 0u 2.95 ±0.4 1.57 3.24 ±0.44 2.32 1.71 

1ZL 



^i°/2 (10 25 y) 5.6 ±1.5 18.7 4.6 ± 1.2 8.54 16.5 



used in the neutrino propagator. The effects of the finite nucleon size and higher-order weak 
currents were taken into account. The two-nucleon short-range correlations (SRC) were 
treated in an extended Brueckner theory (CCM) in a modern self-consistent way, see [25[ 
and Sect. II, that leads to a change in the NMEs M 0v only by a few percents, much less than 
the traditional Jastrow-type representation of the SRC does. 

An important cross-check of the calculations is provided b y a comparison of the present 



results in the spherical limit with the previous ones of Refs. [ll|, [12|, |25j. Though formally 
the adiabatic Bohr-Mottelson approximation is not applicable in the limit of vanishing de- 
formation, it is easy to see that the basic equations fl33l) - fl371) do have the correct spherical 
limit. Details of such a comparison can be found in Ref. |36[, and an excellent agreement 
between the NMEs calculated by the genuine spherical code and the deformed code in the 
spherical limit was found. Also, the partial contributions M^A" 71 ") of different intermediate 
K n states to M 0u for the decay 150 Nd— s- 150 Sm were analyzed in Ref. 3^ . 



The final results for the NMEs, corresponding to the modern self-consistent treatment 
of the SRC O), for Oi//3/3 decays 76 Ge^ 76 Se, 150 Nd^ 150 Sm, 160 Gd^ 160 Dy are listed in 
Table ILTTl As explained in Ref. [36], the difference between the spherical and deformed results 
mainly come from the BSC overlap between the ground states of the initial and final nuclei. 
The strongest effect of deformation on M 0u (the suppression by about 40% as compared 
to our previous QRPA result obtained with neglect of deformation) is found in the case of 
150 Nd. This suppression can be traced back to a rather large difference in deformations of 
the ground states of 150 Nd and 150 Sm. As for the g pp dependence of the Oi//3/3-decay NMEs, 
it is much less pronounced than the dependence of the amplitude of 2v(3f3 decay. A marked 
reduction of the total M' 0l/ for the quenched value of qa can be traced back to a smaller 
pref actor (g A /1.25) 2 in the definition of M' 0u ®. 

In Table HVl the NMEs M 0u for 150 Nd calculated by other approaches are compared. The 
NMEs M 0v for 150 Nd, obtained within the state-of-the-art QRPA approach that accounts for 
nuclear deformation [36j], compares well with the results of the IBM jl6[ and PHFB The 
calculated 0^/9/3-decay half-life Ty 2 corresponding to the Majorana neutrino mass (mpp) = 
50 meV seems to be short enough to hope that the SNO± experiment will be able to approach 
the inverse hierarchy of the neutrino mass spectrum. 



IV. ON THE POSSIBILITY TO MEASURE FERMI Oz//3/3-DECAY NUCLEAR 
MATRIX ELEMENT 

Although there has been great progress in the calculations of the NME M 0v over the last 
decade, but still there can be a substantial scatter in the calculated M 0u by different groups. 
Even more striking, up to a factor of 5, can be the difference in the Fermi part Mp u of the 
total M 0u . 
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Therefore, it would be very important to find a possibility to determine M 0u experimen- 
tally. Partial one-leg transition amplitudes to the intermediate 1 + states have been measured 
by charge-exchange reactions in many nuclei (see |40J and references therein), thereby pro- 
viding important spectroscopic information. However, an attempt to reconstruct the nuclear 
amplitude M 2u of two-neutrino /3/3 decay from the measured amplitudes suffers from large 
inherent uncertainties stemming from unknown relative phases of different intermediate- 
state contributions. Thus, only if a transition via a single intermediate 1 + state dominates 
M 2u , M 2u can consistently be determined. Trying the same way to reconstruct M 0l/ seems 
even more hopeless, since many intermediate states of different multipolarities (with a rather 
moderate contribution of the 1 + states) are virtually populated in the 0^(3(3 decay due to 
a large momentum of the exchanged virtual neutrino. In addition, the transition operators 
in a charge-exchange reaction and Ov/3/3 decay become more and more different for higher 
spins of the intermediate states. 

A pro posal suggesting a way of a direct measurement of was put forward in a recent 
work [4l|. It exploits the similarity between the Fermi part of the neutrino potential in 0v(3(3 
decay and the radial dependence of the Coulomb interaction. The latter is well-known to be 
the leading source of the isospin breaking in nuclei 42, 43]. As shown in Ref. (4lJ, the Fermi 



matrix element M F U can be related to the Coulomb mixing matrix element between the ideal 
double isobaric analog state (DIAS) of the ground state (g.s.) of the initial nucleus 1 and 
the g.s. of the final nucleus. As a result of the Coulomb mixing, the single Fermi transition 
matrix element (Of\T~\IAS) between the isobaric analog state (IAS) of the g.s. of the 
initial nucleus and the g.s. of the final nucleus becomes non-zero. Thus, having measured 
this single Fermi transition matrix element (Of\T~\IAS) by charge-exchange reactions, the 
(V/3/3-decay matrix element M F V can be reconstructed. 

Of course, by measuring only M F V one would not get the total matrix element M 0v 
but rather its sub-dominant part contributing about 20-30% to M 0u . However, knowledge 
of M F U itself brings a very important piece of information. For instance, it can help to 
discriminate between different nuclear structure models in which calculated M F V may differ 
by as much as a factor of 5. More importantly, the ratio M F V /Mq^ may be more reliably 
calculable in different models than M F V and M(Y T separately. Simple arguments put forward 
in Ref. j4l| showed that an estimate M < q T /M f u w —2.5 should hold in a realistic calculation 



(QRPA results |ll|, |12|, |25| do agree with this simple estimate). 

The master relation, derived in Ref. [4l[ in the closure approximation 2 by making use of 
the isospin symmetry of strong interaction H s t r , represents the matrix element M F V in the 
form of an energy- weighted double Fermi transition matrix element: 

M°/ = -^^s(0 f \f-\0i)(0i\f-\0 t }. (38) 

e s 

Here, T~ = J2 a T a * s the isospin lowering operator, the sum runs over all + states of the 
intermediate nucleus z+ ^E1jv-i, w s = E s — (E . + E 0f )/2 is the excitation energy of the 
intermediate state s relative to the mean energy of g.s. of the initial and final nucleus. 



1 This ideal DIAS would be an exact nuclear state if the isospin symmetry were exact. 

2 Using closure of the states of the intermediate nucleus z+ ^ELy_i which are virtually excited in /3/3-decay 
would be an exact procedure if there were no energy dependence in the Qv(3j3 transition operator. A 
weak energy dependence of the operator leads in reality to a "beyond-closure" correction to the total M 0v 
which does not exceed 10%. 
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To account for the isospin-breaking part of H s t r , the r.h.s. of Eq. (|38|) should be slightly 
modified. It is well known that the isospin-breaking terms in H str are in fact fairly small, at 



the level of 2%-3% [42j,|43j, and we arrived in Ref. [41] at the conclusion that the contribution 
of this source of the isospin non- conservation to Eq. f[3"5j) is about 20-30 % of that caused 
by the Coulomb interaction. 

As argued in Ref. [4lj|, the expression ( )38|) in the leading order of the Coulomb mixing 
must be dominated by the amplitude of the double Fermi transition from the initial g.s. via 
its IAS into the final g.s.: 

M f ~ ~^u;iAs(0f\f-\IAS)(IAS\f-\0 l ), (39) 

Here, the second Fermi transition amplitude is non-vanishing due to an admixture of the 
ideal double IAS (DIAS) wave function \DIAS) = , {T ~ ? |0+) in the g.s. of the final 

nucleus: (0 f \f-\IAS) = (O f \DIAS)(DIAS\f-\IAS), and T = (N — Z)/2 is the isospin of 
the g.s. of the initial nucleus. 

In Eq. fl39D, the first-leg matrix element (IAS\f-\0 t ) « yfW Q = \/N - Z and the IAS 
energy ujas are very accurately known. Thus, the total Mp u can be reconstructed accord- 
ing to Eq. (I39p . if one is able to measure the AT = 2 isospin-forbidden matrix element 
(/AS'|T + |0/), for instance in charge-exchange reactions of the (ra,p)-type. 

From the value of Mjjf calculated in a model, the magnitude of the matrix element 
(IAS\T + \0f) can be estimated by using a transformed version of Eq. (!39j) : 

<MS|t+|0 '> = (40) 

Using recent QRPA calculation results for Mp" [llj . this matrix element can roughly be 
estimated as (IAS\T + \0f) ~ 0.005, i.e. about thousand times smaller than the first-leg 
matrix element (IAS\T~\0i). This strong suppression of (IAS\T + \0f) reflects the smallness 
of the isospin-breaking effects in nuclei. 

The IAS has been observed as a pronounced and extremely narrow resonance, and 
its various features have well been studied by means of (p,n), ( 3 He,t) and other charge- 
exchange reactions on the g.s. of a mother nucleus. In this case the reaction cross-section 
at the zero scattering angle can be shown to be proportional to a large Fermi matrix el- 



ement (IAS\T |0j) ~ y/N — Z 46]. Extraction of a strongly suppressed matrix element 
(IAS\T + \0f) from a tiny cross-section of the (n,p)-type reactions on the final nucleus might 
only be possible if there exists a similar proportionality in the (n,p) channel. 

Therefore, a detailed realistic analysis of the corresponding reaction mechanism is needed 
to assess the possibility of extraction of the matrix element (IAS\T + \0f) from the corre- 
sponding reaction cross-sections. 

A first preliminary assessment of the (n,p) reaction at the zero scattering angle was 
done in Ref. [3]. The IAS was treated as a single, well- isolated, state as it appears in 
rather light nuclei. In such a case the Coulomb mixing could be treated perturbatively that 
significantly simplified the consideration. In heavier nuclei the spread of the IAS becomes 
rather significant and should be taken into consideration. 

As argued in Ref. the isospin of the projectile should not be larger than T = 1/2. 
Indeed, the main components of the wave functions |0/) and \IAS) have the total isospin 
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different by two units. Therefore, already for a projectile with isospin T = 1 a common 
entrance and exit isospin channel exists arising from a recoupling of the isospin T = 1 of 
the projectile with the main components of the wave functions of the target and daughter 
nuclei. In such a case extraction of the information about small isospin impurities from the 
corresponding reaction cross-section seems barely possible. 

Thus, the only probes which seem feasible are of the isospin T = 1/2 ((n,p), (t, 3 He), 
...reactions). However, it is still not guaranteed that the reaction cross-section cr(0j — >• 

IAS) for these probes is proportional to a strongly suppressed matrix element (IAS\T + \0f) , 
since the other isospin impurities in the wave functions |0/) and \IAS) may have a larger 
effect on the reaction cross-section. 

A preliminary assessment of the (n,p) reaction at the zero scattering angle in the 
aforementioned perturbative mixing approximation was performed in Ref. [44] for an 
intermediate-mass /3/3-decaying nucleus 82 Se. In fact, it was shown that the tiny cross- 
section cr np (0j — > IAS) is indeed dominated by the admixture of the DIAS in the g.s. of the 
final nucleus. However, the spread of the IAS of 82 Se is rather significant that may question 
the perturbative treatment of the isospin mixing. 

A /3/3-decaying nucleus in which such a perturbative treatment may be justified in reality 
is 48 Ca. This case was considered in detail in Ref. (45| . 

The IAS of 48 Ca is a state with J 71 = + ,T = 4, T z = 3 at the excitation energy of 
E x =6.678 MeV {ojias ~8.5 MeV) in 48 Sc. It lies under the threshold of particle emission 
and with almost 100% probability decays to 1 + state at E x =2.517 MeV by the emission of a 
7-quantum with i? 7 =4.160 MeV j47|. This 7-decay energy is much higher than the 7-decay 
energies from + states with normal isospin T = T z = 3 surrounding the IAS (which decay 
by a cascade), and could be used as a unique experimental tag telling that the IAS indeed 
was excited in a reaction. 

There are strong arguments that the IAS of 48 Ca must be a single state without frag- 
mentation. The state-of-the-art measurement of 48 Ca( 3 He,t) 48 Sc(IAS) reaction 48j does in 
fact contribute to clarification of this issue as discussed below. 

Fragmentation of the IAS may occur only if there are several + states with the normal 
isospin around the IAS to which the IAS may strongly couple. In other words, the total 
number of the + states within the IAS spreading width must be greater than one (for 
nuclei around A = 50 T\ is typically about few keV). Based on the back-shifted Fermi-gas 



model 49j, the mean level spacing between the + states of the normal isospin in the vicinity 



of the IAS in 48 Sc was estimated in Ref. |45j to be about 50-70 keV. Then, if the IAS were 



essentially spread over those + states, the experiment [48J would have been able to resolve 

components of the IAS fine structure. The fact that no fine structure was observed can 

easily be understood from a comparison of a typical of the order of few keV with a much 
larger mean level spacing. 



The cross-section for the reaction 48 Ti(n,p) 48 Sc(IAS) was estimated in Ref. 45J to be 



4^]§ ~ 20 nb/ (sr MeV) for the energy of the incident neutrons around 100 MeV. Note, that 
by choosing a smaller neutron incident energy this estimate can further be improved by a 
factor of 2-3, due to the increasing Fermi unit cross-section |46] . The Coulomb mixing of 
the IAS with the isovector monopole resonance was estimated to modify the above value of 
by few percents, and therefore can be neglected. 
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V. THE RESONANT NEUTRINOLESS DOUBLE ELECTRON CAPTURE 



The resonant Oz/ECEC (neutrinoless double electron capture) was considered as process, 
which might prove the Majorana nature of neutrinos and the violation of the total lepton 
number, by Winter [50| already in 1955. The possibility of a resonant enhancement of the 
Oz/ECEC in case of a mass degeneracy between the initial and final atoms was pointed out 



by Bernabeu, De Rujula, and Jarlskog as well as by Vergados about 30 years ago [51|, |52 
They estimated the half-life of the process by introducing different simplifications: i) non- 
relativistic atomic wave functions at nuclear origin; ii) qualitative evaluation of NME of 
the process; iii) the degeneracy parameter A = Ma,z — M A * Z _ 2 was assumed to be within 
the range (0,10) keV representing the accuracy of atomic mass measurement at that time. 
Ma,z and M A * Z _ 2 are masses of the initial and final excited atoms, respectively. A list of 
promising isotopes based on the degeneracy requirement associated with arbitrary nuclear 
excitation and on the natural abundance of daughter atom was presented. 

In 2004 Sujkowski and Wycech [53| and Lukaszuk et al. (54j analyzed the Oz/ECEC pro- 
cess for nuclear + — > + transitions accompanied by a photon emission in the resonance and 
non-resonance modes. By assuming \mpp\ = 1 eV and 1 a error in the atomic mass deter- 
mination the resonant Oz/ECEC rates of six selected isotopes were calculated by considering 
the perturbation theory approach. 

In 2009 a new theoretical approach to the Oz/ECEC, a unified description of the oscilla- 
tions of stable and quasistationary atoms, was developed by Simkovic and Krivoruchenko 
0, HEJ. A comprehensive theoretical study of this process for the light Majorana neutrino 
mass mechanism was performed jsj. It was shown that effects associated with the relativistic 
structure of the electron shells reduce the Oz/ECEC half-lives by almost one order of magni- 
tude and that the capture of electrons from the npi/ 2 states is only moderately suppressed 
in comparison with the capture from the ns\/2 states unlike in the non-relativistic theory. 
Selection rules for associated nuclear transitions were presented saying that a change in the 
nuclear spin J > 2 are strongly suppressed. New transitions due to the violation of parity 
in the Oz/ECEC process were proposed, e.g., nuclear transitions + — > , 1 are compatible 
with a mixed capture of s- and p-wave electrons. Based on the most recent data and realistic 
evaluation of the decay half-lives, a complete list of the most perspective isotopes for which 
the Oz/ECEC capture may have the resonance enhancement was provided, for further exper- 
imental study. It inludes 96 Ru, 106 Cd, 124 Xe, 136 Ce, 152 Gd, 156 Dy, im Er, 168 Yb, im W, 184 0s 
and 190 Pt 0]. By assuming \m,pp\ = 50 meV and an appropriate value of NME, half-lives of 
some of the isotopes were found to be as low as 10 25 years in the unitary limit. It is about 
one order of magnitude shorter than the Oz//3/3 half-life of 76 Ge for the same value of the 
effective mass of Majorana neutrinos. 

The inverse value of the half-life of resonant neutrinoless double electron capture 



In 2 \V ab (r)\ 2 



T ab (41) 



T /f CEC (j7r) A 2 + l T 2 b 

where J n denotes angular momentum and parity of final nucleus. The degeneracy parameter 
can be expressed as 

A = M AZ - M* A * Z -2 = Q~B ab - £ 7 , (42) 

where Q stands for a difference between the initial and final atomic masses in ground states 
and E 1 is an excitation energy of the daughter nucleus. B ab = E a + E b + Eq is the energy 
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of two electron holes, whose quantum numbers (n, j, I) are denoted by indices a and b and 
Ec is the interaction energy of the two holes. The binding energies of single electron holes 



E a are known with accuracy with few eV |56j. The width of the excited final atom with the 
electron holes is given by 

r ab = r a + r b + r*. (43) 

Here, T a ^ is one-hole atomic width and T* is the de-excitation width of daughter nucleus, 
which can be neglected. Numerical values of T a b are about up to few tens eV. By factor- 
izing the electron shell structure and nuclear matrix element for lepton number violating 
amplitude associated with nuclear transitions + — > J w = ±:L , 1 ±:L one gets 



V ab (r) = — Glm^-^L < F ab > M^^'(r). (44) 



OuECECi 



Here, < F a b > is a combination of averaged upper and lower bispinor components of the 
atomic electron wave functions i and M^ ECEC {r) is the nuclear matrix element. We note 
that by neglecting the lower bispinor components M 0l/ECEC (0 + ) takes the form of the 0v(3(3- 
decay NME for ground state to ground state transition after replacing isospin operators r~ 
by r + . R is the nuclear radius and is the axial- vector coupling constant. 

The probability of the Oz/ECEC is increased by many orders of magnitude provided the 
resonance condition is satisfied within a few tens of electron-volts. For a long time there was 
no way to identify promising isotopes for experimental search of Oz/ECEC, because of poor 
experimental accuracy of measurement of Q-values of the order of 1 - 10 keV for medium 



heavy nuclei. Progress in precision measurement of atomic masses with Penning traps [57 



59| has revived the interest in the old idea on the resonance OzaECEC. The accuracy of 



Q-values at around 100 eV was achieved. The estimates of the OzaECEC half-lives were 



recently improved by more accurate measurements of Q- values for 74 Se [60, |6jj , 96 Ru [6£_ 
106 Cd @,pl, 102 Pd B, 112 Sn H, 120 Te [^i^Ce @, 144 Sm Q, 152 Gd 0, 156 Dy 
@, 162 Er 19, 164 Er Q, 168 Yb @ and 180 W [H|. It allowed to exclude some of isotopes 
from the list of the most promising candidates (e.g., 112 Sn, 164 Er, 180 WW) for searching the 
Oz/ECEC. 

Among the promising isotopes, 152 Gd has likely resonance transitions to the + ground 
states of the final nucleus as it follows from improved measurement of Q-value for this 
transition with accuracy of about 100 eV [671. A detailed calculation of the Oz/ECEC of 152 Gd 



was performed in [7_2j (see Table IHlj) . The atomic electron wave functions were treated in the 
relativistic Dirac-Hartree-Fock approximation The NME for ground state to ground 

state transition 152 Gd — > 152 Sm was calculated within the proton- neutron deformed QRPA 



with a realistic residual interaction [67[. For the favored capture of electrons from K and 
L shells in the case of 152 Gd the 0i/ECEC half-life is in the range 4.7 x 10 28 - 4.8 x 10 29 
years. This transition is still rather far from the resonant level. Currently, the Oz/ECEC 
half-life of 152 Gd is 2-3 orders of magnitude longer than the half-life of decay of 76 Ge 

corresponding to the same value of \mpp\ and is the smallest known half-life among known 
Oz/ECEC. 

The resonant Oz/ECEC has some important advantages with respect to experimental 
signatures and background conditions. The ground state to ground state resonant Oz/ECEC 
transitions can be detected by monitoring the X rays or Auger electrons emitted from excited 
electron shell of the atom. This can be achieved, e.g., by calorimetric measurements. The 
de-excitation of the final excited nucleus proceeds in most cases through a cascade of easy 
to detect rays. A coincidence setup can cut down any background rate right from the 
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beginning, thereby requiring significantly less active or passive shielding. A clear detection 
of these 7 rays would already signal the resonant Oz/ECEC without any doubt, as there 
are no background processes feeding those particular nuclear levels. We note that standard 
model allowed double electron capture with emission of two neutrinos, 

e 6 - + e 6 - + (A, Z) ->(A,Z- 2)** + z/ e + v e , (45) 



is strongly suppressed due to almost vanishing phase space [52 . 

Till now, the most stringent limit on the resonant Oz/ECEC were established for 74 Se 
Cd 1 76[ and 112 Sn (7?]]. However, following recent theoretical analysis [3] none of 



these resonant Oz/ECEC transitions is favored in the case of light neutrino mass mechanism. 
The ground state of 74 Se is almost degenerate to the second excited state at 1204 keV in 



the daughter nucleus 74 Ge, which is a 2 + state [78| and is disfavored by the selection rule 
[3|. The TGV experiment situated in Modane established limit on the Oz/ECEC half-life of 
1.1 x 10 20 years [76J. Subject of interest was the Oz/ECEC resonant decay mode of 106 Cd 
(KL-capture) to the excited 2741 keV state of 106 Pd. For a long time the spin value of this 
final state was unknown and it was assumed to be J 71 " = (1,2) + . However, a new value for 
the spin of the 2741 keV level in 106 Pd is J = 4 + and this transition is disfavored again 
due to selection rule. A search for the resonant Oz/ECEC in 106 Cd was carried out also at 
the Gran Sasso National Laboratories with the help of a 106 CdWO^ crystal scintillator (215 
g) enriched in 106 Cd up to 66% @. It was found that the resonant Oz/ECEC to the 2718 
keV (J* is unknown), 2741 keV (J* = 4+) and 2748 (J* = (2,3)") keV excited states of 
106 Pd are restricted to T^ ECEC > 4.3 x 10 20 yr (KK-capture), T^ ECEC > 9.5 x 10 20 yr 
(KLi-capture) and Ty EOEC > 4.3 x 10 20 yr (KL 3 -capture), respectively. We note that the 
2718 excited state 7 decays by 100% into the 3+ state at 1557.68 keV state, which again 
excludes a possibility of J = 0, 1 for this state. Further, we already mentioned above that 
a new mass measurement [(mJ has excluded a complete mass degeneracy for a 112 Sn decay 
and has therefore disfavored significant resonant enhancement of the Oz/ECEC mode for 
this transition. Recently, a first bound on the resonant Oz/ECEC half-life of 136 Ce of about 
1.1 x 10 15 years was measured [80] . 



VI. SUMMARY 

Many new projects for measurements of the Oz//3/9-decay have been proposed, which hope 
to probe effective neutrino mass down to 10-50 meV. An uncontroversial detection of 
the Oz//3/3decay will prove the total lepton number to be broken in nature, and neutrinos to 
be Majorana particles. There is a general consensus that a measurement of the Oz//3/3-decay 
in one isotope does not allow us to determine the underlying physics mechanism. It is very 
desired that experiments involving as many different targets as possible to be pursued. There 
is also a revived interest to theoretical and experimental study of the resonant Oz/ECEC 
(neutrinoless double electron capture), which can probe the Majorana nature of neutrinos 
and the neutrino mass scale as well. The Oz/ECEC half-lives might be comparable to the 
shortest half-lives of the Oz//3/3 decays of nuclei provided the resonance condition is matched 
with an accuracy of tens of electron-volts. There is a lot of theoretical and experiment effort 
to determine the best Oz/ECEC candidate. 

Nuclear matrix elements of these two lepton number violating processes need to be eval- 
uated with uncertainty of less than 30% to establish the neutrino mass spectrum and CP 
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violating phases. Recently, there has been significant progress in understanding the source 
of the spread of calculated NMEs. Nevertheless, there is no consensus among nuclear the- 
orists about their correct values, and corresponding uncertainty. The improvement of the 
calculation of the nuclear matrix elements is a very important and challenging problem. We 
presented improved calculation of the Oz//3/3-decay and Oz/ECEC NMEs, which includes a 
consistent treatment of the two-nucleon short-range correlations and deformation effects. In 
addition, a possibility to measure the 0^/9/9-decay NME was addressed. 
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